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Abstract 

The stability of a /3-sheeted conformation and its transition into a random coil are studied with 
a 2D lattice biopolymer model. At low temperature and low external force, the polymer folds back 
and forth on itself and forms a /3-sheet. Our analytical calculation and Monte Carlo simulation 
reveal that a co-operative /3-sheet-random coil transition takes places when the temperature or 
force is increased, with a dramatic decrease in the contact number. These predictions are in 
good agreement with experiments on titin protein. This transition is not a real phase-transition, 
indicating that backbone hydrogen-bonding alone is unable to stabilize a distinct /3-sheet phase. 
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Single-molecule manipulation techniques were used by Gaub and coworkers and several 
other groups [|1], to study the unfolding of proteins in real time. A co-operative "all-or- 
nothing" transition between globular and randomly coiled states was observed at a pulling 
force of ~ 10 2 piconewtons. A stretched coil segment can be modelled as a wormlike chain of 
persistence length ~ 10 _1 nanometers [^J, [J, and several (atomic-level) molecular dynamics 
(MD) simulations have studied the initial stage of unfolding and the final stage of refolding 
(see, e.g., || |J). Monte Carlo (MC) simulations || of lattice models also have been used to 
study the unfolding kinetics. The protein used in the experiments, namely titin imunoglob- 
ulin domains, is rich in /3-sheets Jl|, 0], and numerical work demonstrated that the transition 
to a random coil was initiated by the rupture of only a few hydrogen bonds between adjacent 
/3-strands 0. Theoretically, however, it is still unclear why the globule-coil transition is a 
highly co-operative process, and why there are dramatic differences between the unfolding 
behavior of helix-dominated and /3-sheet-dominated proteins || . 

An issue of current interest is to study the generic statistical mechanics of biopolymer (3- 
sheet states. In this paper, the deformation and unfolding of /3-sheets are addressed within 
an exactly-solvable 2D model. Because of the attractive interaction between monomers 
along the chain, the polymer could fold back and forth on itself, thereby maximizing its 
contact energy. However, a /3-sheet, which is a relatively ordered structure and is favorable 
in energy, is deficit in configurational entropy compared with an irregular coil configuration. 
Because of this competition between energy and entropy, it is demonstrated by our analytical 
and simulational calculations that, when the temperature or the external force approaches 
some critical value, a small variation in temperature or force could lead to a dramatic 
change in the polymer's microscopic configuration. This behavior corresponds well both 
with experimental observations and MD simulations on titin. However, we suggest that this 
highly co-operative transition is not a phase-transition in the strict sense; the nonspecific 
hydrogen-bonding interactions between the backbone molecules of a protein alone are hence 
unlikely to be able to stabilize a distinct /5-sheet phase. 

The mathematical model is constructed as follows: Consider a biopolymer of N identical 
units (monomers) located on a square lattice (Fig. [[]). The distance between two consecutive 
monomers is do, and to facilitate the calculation the chain is assumed to be directed, namely 
the bond between monomers i and i + 1 can not point to the — z direction (Fig. [I]). If 
two monomers with indices i and i + m (m > 3) occupy nearest neighbor lattice sites, an 
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attractive energy of strength e is gained and these two monomers are then "in contact". 
Each monomer can be in one of three states, with respectively zero, one, or two contacts. 
By folding back and forth on itself, the polymer can gain contact energy, and a segment 
of m columns in which there is at least one pair of monomers in contact between any two 
adjacent columns is referred to as a "/3-sheet". Two /5-sheet segments are separated by 
a "coil segment" of m columns (m > 0), in which all the monomers are free of contacts 
(Fig. |). 

We first discuss the free energy of the /3-sheet state. Under the action of a force / along 
direction z , the energy of a /3-sheet of m columns is — e Yh=^ v (^> h+i) — m f a o, where 
v(/j, l i+ i) = Li — 1 (if li < and Zj + i — 1 otherwise; /j is the number of monomers in the 
i-th column. The partition function (PF) Zp{n) of a /3-sheet of length n is 

[n/2] m -l 

Mn) = E E • • • E J m II « v(Wl+l) , (i) 

m=2Zi=2 l m =2 i=l 

where 51 — 1 (i — j) and (otherwise); a = exp(/3e), b = exp(/3/a ), and (3 = 1/T, T being 
the temperature (the Boltzmann constant is set to 1). 

There can be many-body interactions among the monomers in a /3-sheet and a contact 
pair can be formed between two remote monomers along the chain, consequently, direct 
evaluation of Zp(ri) for a large system is difficult. To circumvent this, we define an auxiliary 
function T m ((,k) = J2 h=2 ■ ■ ■ £« m=2 ( k ^a v ^M~ h \[™S 2 l (C / 'a) u ^ k ' k + 1 ' \( , 'of™ for m > 2 and 
Ti(C, k) = ( 2k /a{a - ( k ). One finds 

C^UCl) M C k (l-a)T m -i(C,k + l) 
mlC " } ~ a(l-C k ) a(a -C k )(l-C k ) ' K ) 

Define Lp^, k) = J2m=i ^ m ^m(C) k), and from the recursive relation Eq. (fj) one obtains the 

following two equations: (1) L P (C,2) = A{() + B(C)L P (C,1) and (2) L p (C,2) = H{() + 

J(C)L /3 (^, 1). After eliminating Lg(£,2) from these two equations, one finds that 

MC;1) - j(o-B(cy (3) 

In Eq. ©, A(C) = C(l-C)/(o-l), B(C) = ( a -()(be + a(-a)/b(a-l)(, H = {b/a)[CKa- 

C 2 ) + E J =3nS0(C,jX 2 V(a-CO],andJ=(Va)[CV(l-C 2 )+E^ 3 n}= 1 2 0(CjX 2 V(l^ 
where 0(C, j) = 6(1 - a)C, k /a(a - C fc )(l - C fc )- 

Both H(Q and J(() are convergent when ( e [Ojl]- When A = 1 — ( — ► + , our 
analytical and numerical calculations suggest that A(Q — H(C) ~ > exp(— b/aA)f 1 (A) and 



J(Q - B(C) - exp(-6/aA)/ 2 (A), with /i(A)// a (A) - 6/a(a - 1)- Here /i(A) and / 2 (A) 
are two polynomials of A which contain terms with negative powers of A. 

The grand partition function (GPF) of the /3-sheet, Gg(C) = J2n=4(( / a ) n Z p(n) , is nothing 
but Lp((, 1) — b 2 ( 2 /a(a — (). Gp(Q could be calculated based on Eq. (^) at given (. It keeps 
increasing as ( approaches (p{< 1) from and is divergent at (p, where (p is the smallest 
positive root of the equation J(() — B(() = (see Fig. The free energy density of 

a /5-sheet segment equals to pp = fc^Tln^/a). The divergence of Gp at (p indicates the 
absence of a distinct /3-sheet phase. 

Armed with the partition function of a /5-sheet, we proceed to investigate its stability 
and the configurational transition of the whole polymer system. To this end, one needs to 
know the free energy of a coil segment 0. The total PF of a coil segment of length n under 
the action of an external force (/, incorporated in the parameter b) is expressed as 

n n 

ZUn) = 2 £ ZUn, m)b m + £ Z^(n, m)b m (4) 

m=l m=l 

and Z coi i(0) = 1. In this expression, Z^ oil (n, m) is by definition the total number of allowed 
configurations for a coil of length n whose two ends are separated by m lattice steps along 
z and whose last column has only a single monomer; Z™ n (n, m) is defined similarly, but 
with the last column containing at least two monomers. The coefficient 2 in Eq. (Q) appears 
because the first column of a /3-sheet following a coil segment that ends with a column 
containing a single monomer can have two possible orientations (Fig. |l|). 

One finds Z^ oil (n, m) = Zl Q - lX (n— 1, m— l)+Z™ il (n— 1, m— 1), and Z™ i{ (n, m) = Z^ oil (n,m+ 
1) + Z s coil (n - 2, m - 1), with Z c s oil (0, m) = 5° m , Z™ a (0, m) = Z^{1, m) = 0, 2^(2, m) = 5 x m , 
and Z™ a (n, 1) = 1 — 5^. With these preparations the GPF of the coiled state is calculated 
to be 

r rn _v r n-n 7 ,s_ a(a-()(a + b() 

( ( )_ a 3 -a2(6 + l)C-&V> (5) 

It follows that the free energy density of a coil segment equals p coi i = kBT\n(( coi i/a). Ccoii 
is the point where G coi i(C) = +oo: Ccoii = ln[(a/6)(( v / A + 6/2) 1 / 3 - (V\ - 6/2) 1 / 3 )] with 
A = (6 + l) 3 /27 + 6 2 /4- 

In the lattice model, any configuration of the polymer is a chain of /3-sheets and random 
coils, with the two conformations occurring alternately (Fig. [I]). The total PF is then 

[N/4] 

z(n) = EEEE-EEE 

s=0 i =0ii=4ii=0 j s =4i s =Ojo=0 
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fc=0 

Finally the GPF of the whole system may be calculated: 

CrtotalKJ - 2^ C O Ztotal^VJ - — -. (7 J 

Ar=0 1 ~~ ^ lU^coil [ U 

Because of the facts that Gp((,p) = +00 and G C oii(Ccoii) = +00, it is always possible to 
find a root (Ctotai) of the equation 

G'coil(Ctotal)G r / 3(C t otal) = 1 (8) 

in the range of < Ctotai < ixiin(^" /3 , Ccoii)- Consequently, the free energy density of the whole 
system is neither that of the /3-sheet state nor that of the coil state: ptotai = ksT ln(Ctotai/o)- 
There is no singularity in the total free energy expression, and hence no real phase-transition 
process. Almost forty years ago, it was revealed by intensified theoretical investigations that 
the protein a-helix-coil transition is not a phase-transition, provided that only the main- 
chain hydrogen-bonding interactions are considered (see, e.g., || [|); we are surprised to find 
that the same conclusion holds for the /3-sheet transition. Although there is no strict phase- 
transition between /3-sheet and random coil in the model, our analytical and simulational 
results presented below indicate that, when temperature and the external force are low 
enough the polymer is essentially in the form of a /3-sheet, and there is an rapid change in 
the polymer's total extension and its total number of contacted pairs when the temperature 
or the force is increased. 

The total contact number a (scaled by N) could be calculated by a = 1 — d ln((totai) /din a 
and the total extension £ (scaled by Nao) could be calculated by £ = —bd\ia((totai)/db. 
When there is no external force (/ = 0), the temperature-contact number profile and the 
temperature-extension profile are shown in Fig. [5| And in Fig. |] the relationship between 
force and extension and between force and contact number are shown at a fixed temperature 
of T = 0.591e. When the temperature and the force is small, Ctotai approaches (the 
divergent radius of Gp) and its value is difficult to be numerically determined in much 
precision through Eq. @. Therefore, the analytical calculation in Fig. |3| is shown only for 
T > 0.8e and that in Fig. f| for / > 0.6e/ao- In Figs. and f| the data obtained by a MC 
simulation is also shown. The Monte Carlo simulation is performed with a polymer size 
of N = 1024 (the results for iV = 5000 are essentially the same as the data shown here). 
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There are three kinds of elementary updates of configurations: (i) one monomer is moved 
from one column to another column or stacked on the right-hand side of the last column; 
(ii) the positions of two columns are exchanged; (iii) the direction of one column is changed 
by 180°. Special attention is paid to make sure that ergodicity and detailed balance are 
obeyed. Each simulational data is obtained through 10 9 MC steps (it takes about 35 hours 
in a Compaq station). The energetically most favorable state is a /3-sheet whose width and 
height both equal to y/N . Therefore, in the numerical simulation data shown in Figs. ||] and 
H for a chain of N = 1024, one has o < 0.935 and £ > 0.0313. For an infinite system, a — ► 1 
and t — > as T — > 0. 

Figure [| indicates that when / = and T < 0.75e almost all the monomers of the polymer 
are in the /3-sheet state. At T ~ 0.75e there is a plateau in the T-o and T—l profiles, 
suggesting that a co-operative transition into random coiled structure takes place. During 
this process, the monomers lose contacts and the polymer is extended. As temperature is 
further increased, the configuration of the polymer becomes essentially random and there are 
only small changes of the total extension and the total contact number. Similar co-operative 
/3-sheet-coil transition is observed when the temperature is fixed but the external force keeps 
increasing. At T = 0.591e Fig. ^ indicates the transition takes place at / ~ 0.5e/cto. In the 
experimental investigations of Gaub and coworkers [|T| a co-operative globule-extended coil 
transition is observed for the /3-sheet dominated protein titin; and in a later experiment of 
Fernandez and coworkers it was suggested that the force-induced denaturation of protein 
titin should follow the same paths as the temperature-induced denaturation of this protein. 
The present theoretical work seems to confirm both of these experimental observations: 
a co-operative transition from /3-sheet to extended coil will occur either by increase the 
environmental temperature or by increase the mechanical force acting on the two ends of 
the polymer. 

Conformational transitions in single biopolymers have attracted the interest of theorists 
for quite a long time, and the introduction of simple mathematical models have brought a 
deep understanding of the issues of the protein a-helix-to-coil transition ||, double-stranded 
DNA melting JTIJ, /3-hairpin formation in oligopeptides |IT|, and RNA secondary structure 



denaturation by heat pi or by force [I3]|, etc.. The current work extends these efforts to a 
new model system, namely, /3-sheet formation in biopolymers. It suggests that a co-operative 
/3-sheet-coil transition will happen at certain temperature or mechanical stress. To make it 
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more realistic, other effects such as the sequence heterogeneity and the interactions between 
adjacent layers of /^-sheets || could be considered. On the other hand, the results presented 
here show that even without these complications the structural transition of a /3-sheet is 
"two-state-like" . 

Recently, the MC work in [H suggested that pulling a collapsed polymer in 2D leads to a 
co-operative transition into an extended random coil, and the MC work in [[T5j suggested this 
process actually may be a continuous transition. The present effort is also complementary 
to these MC simulations. 

The financial support of the A. v. Humboldt Foundation is acknowledged. The author 
is grateful to Professor R. Lipowsky for his hospitality, and to Dr. Julian Shillcock for a 
careful reading of the manuscript and many good suggestions. 
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FIG. 1: A model directed biopolymer. Black squares are monomers and the arrow denotes the zq 
direction. The configuration shown here has two /3-sheets and two coils. 
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FIG. 2: The value of the partition function Gp(Q as its divergent radius £/? is approached. The 
temperature is set to T = 0.590928e and the force is / = 0.73866e/ao or / = (the inset). 
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FIG. 3: Dependence of the contact number and the extension on the temperature at zero force. 
Symbols are MC simulation data performed for N = 1024. 
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FIG. 4: Dependence of the contact number and the extension on the external force at fixed 
temperature T = 0.590928e. Symbols are MC simulation data performed for N = 1024. 
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